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ABSTRACT 

Aims. We analyse 3-min oscillations of microwave and EUV emission generated at different heights of a sunspot atmosphere, 
studying the amplitude and frequency modulation of the oscillations, and its relationship with the variation of the spatial structure of 
the oscillations. 

Methods. High-resolution data obtained with the Nobeyama Radioheliograph, TRACE and SDO/AIA are analysed with the use of 
the Pixelised Wavelet Filtering and wavelet skeleton techniques. 

Results. 3-min oscillations in sunspots appear in the form of repetitive trains of the duration 8-20 min (13 min in average). The 
typical interval between the trains is 30-50 min. The oscillation trains are transient in frequency and power We detected a repetitive 
frequency drifts of 3-min oscillations during the development of individual trains. Wavelet analysis shows that there are three types 
of the frequency drift: positive drifts to high frequencies, negative drifts and fluctuations without drift. Negative drifts were found to 
occur more often. The start and end of the drifts coincide with the start time and end of the train. The comparative study of 3-min 
oscillations in the sequences of microwave and EUV images show the appearance of new sources of the oscillations in sunspots 
during the development of the trains. These structures can be interpreted as waveguides that channel upward propagating waves, 
responsible for 3-min oscillations. A possible explanation of the observed properties is the operation of two simultaneous factors: 
dispersive evolution of the upwardly-propagating wave pulses and the non-uniformity of the distribution of the oscillation power over 
the sunspot umbra with different wave sources corresponding to different magnetic flux tubes with different physical conditions and 
line-of-sight angles. 

Key words, sunspot - frequency drifts - EUV emission - oscillation 



1. Introduction 

Analysis of multi-wavelength data obtained with the new gen- 
eration of solar observational tools allows us to study dynamics 
of the solar atmosphere at different heights with high spatial, 
spectral and temporal resolution. One of the interesting topics 
is 3-min oscillations in sunspots (e.g. Orrall [T9661 Elliott [T9691 
Liu 11974) . Fluctuations in the 3-min range are relatively weak 
at the photospheric level of sunspots, but they are significantly 
enhanced in their chromospheres (Lites[1986' Bard and Carlsson 
[19971 Tziotziou et al. 2002, Kobano v et al . 2008 20111 and in 
the transition region (Gurman et al. 119821 Lites 119921 Bogdan 
12000b . Observations of the photospheric magnetic field oscilla- 
tions (Nort on et al. IT9991 Balthasar et al. [19981 Bellot Rubio 
et al. I2OOOI 1 show a fine fragmentation of the magnetic field in 
the regions where the main power of fluctuations is localised 
inside sunspots. For example, optical observations show that 5- 
min oscillations of the magnetic field are mainly concentrated 
in the isolated tubes of magnetic field (pores) outside the big 
sunspots, as well as at the umbra-penumbra boundary (Zhugzhda 
et al. 200()| . Blinkers or explosive events are also observed above 
sunspots with a similar time scale, from 400 to 1600 s (Harrison 
et al. [19991 Brkovich et al. [2000l i. 



The first observations of the fine spatial structure of oscil- 
lating radio sources above strong sunspots at the transition re- 
gion level were obtained by Gelfreikh et al. ( 11999b with the 
Nobeyama Radioheliograph. That study showed that the main 
source of broadband oscillations is concentrated in the centre of 
the source of the polarised radio emission at 17 GHz, coincident 
with the sunspots umbra. Further investigations in the radio band 
(Shibasaki et al., 2001) using the correlation curves and images 
of radio sources showed that 3-min oscillations are present in 
sunspots during long time. It was concluded that the brightness 
fluctuations of the radio sources could be caused by the den- 
sity and temperature perturbations in acoustic waves propagat- 
ing through the third gyroresonant level (corresponding to the 
magnetic field of about 2 kG). 

Nindos et al. (12002b used the difference radio maps obtained 
by VLA at 8.5 and 5 GHz with high spatial resolution. They 
showed that the amplitude of 3-min oscillations varies irregu- 
larly in time, with a rapid rise phase, and a rather slow decay 
phase. The spatial sources of the amplitude variations are lo- 
calised in fine patches of a circular shape. 

Sych and Nakariakov l2008l studied the spatial, temporal and 
phase behaviour of sunspot oscillations, applying the method of 
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Fig. 1. Left panel: Full-disk microwave polarisation image of the Sun at 17 GHz on May 03, 2005 at 02:45 UT. The rectangle shows 
the radio source associated with the sunspot. Right panel: Zoomed image of the radio source shown in the left panel; the dashed 
contours show the L-polarisation image overlapping with the SOHO/MDI magnetogram taken at 01 :35 UT. The grey scale indicates 
the polarity of the magnetic field. 



Pixelised Wavelet Filtering (PWF). It was shown that the source 
of 3-min oscillations is located in the centre of the microwave 
source coinciding with the sunspot umbra. On the other hand, 5- 
min oscillations were found to be mainly located in small, sym- 
metric patches at the the umbra-penumbra boundary, similar to 
the findings obtained with VLA. Oscillations in diff'erent patches 
were found to have different phases. 

Phenomenological relationship between 3-min oscillations 
of the radio flux from a sunspot and quasi-periodic bursts of en- 
ergy release in coronal active regions nearby was established by 
Sych etal_p009 and Nakariakov et al. (20 lO). Recently, Sych 
et al. (12010b used the system for automatic detection of solar 
oscillations (http://pwf.iszf.irk.ru), that showed that the power 
of 5-min oscillations in sunspots is concentrated in ring-shaped 
sources at the upper chromosphere or transition region, which 
coincide with the footpoints of coronal EUV fans. These spatial 
details are associated with both the standing waves (symmetric 
spatial patches) and propagating waves (seen as transient, about 
10-30 min, v-shaped increases in the spatial distribution of the 
oscillation power). The amplitude increase was seen simultane- 
ously in microwaves and in EUV. 

Rendtel et al. (I2003I I detected variable changes (drifts) of the 
periods of sunspot oscillations, using the SOHO/SUMER spec- 
trograph (CIV, NeVIII) and FIT, 171 A . The fluctuations of the 
Doppler shift and the emission intensity occur in the same fre- 
quency range, 3-7 mHz, with a significant time-drift to higher 
frequencies. Most of the fluctuations in the transition region and 
lower corona are concentrated near the frequency of 6.2 mHz, 
and form oscillation trains lasting from 10 to 20 min, some- 
times up to 40 min. A possible explanation of the drifts is the 
assumption that the curvature of the magnetic tubes above the 
sunspot leads to different directions of MHD wave propaga- 
tion relative to the observer and, accordingly, to the periodic 
changes of the Doppler velocity and the sign of the frequency 



drift. Fludra ( I200II) found a significant frequency drift of 5.8 
mHz oscillations above sunspots in the transition region lines 
with SOHO/CDS. The typical duration of the train was found to 
be about 15 minutes. Different frequencies in the spectrum were 
found to have different spatially separated sources. 

Two main mechanisms were suggested to explain 3-min 
sunspot oscillations: a chromospheric cavity (or resonator) and 
oscillations at the acoustic cut-off frequency. The first mech- 
anism implies that the chromosphere is a resonant cavity for 
waves, which are reflected from the photosphere due to sharp 
density gradient and from the transition region due to the sharp 
temperature gradient (Leibacher and Stein II981I Leibacher 
et al. 11982b . It was recently shown numerically that quasi- 
monochromatic acoustic oscillations in the 3-min band are a 
natural response of the chromospheric resonator to a broad- 
band (e.g., impulsive) excitation (201 lb . The second mechanism 
implies that the 3-min oscillations are resulted from the wake 
behind propagating disturbances in the stratified atmosphere, 
which oscillates at the acoustic cut-off frequency (Fleck and 
Schm itz [T99n Kal kofen e t al. 1 19941 Sutmann and Ulmschneider 
119951 Kuridze et al. 12009b . The nonUnear wake behind propagat- 
ing acoustic pulses may develop in quasi-periodic shocks, which 
may shape the plasma properties in the chromosphere and the 
lower corona. These quasi-periodic shocks have been recently 
suggested as a mechanism for the quasi-periodic appearance of 
solar spicules and 5-min oscillations in intensity observed in the 
solar corona (Murawski and Zaq arashviU 120 1 Ol Zaqarashvih et 
al. l20m i. 

In this paper we present results of a comprehensive study 
of 3-min oscillations at different heights over sunspots with the 
use of high-resolution EUV and microwave data obtained with 
SDO/AIA, TRACE and NoRH. The aim of the paper is to reveal 
the fine spatial-temporal structure of 3-min oscillations, con- 
tributing to the understanding of physical mechanisms for their 
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Fig. 2. A correlation curve of NoRH in the circular polarisation 
(R-L) channel. The zoomed time interval, selected by the rectan- pjg, 3, Qjobal wavelet spectrum of the NoRH correlation curve 
gle, shows the smoothed 3-min radio oscillation profile between on May 3, 2005. Arrows indicate the main peaks in the spectrum. 
23:55 and 03:15 UT. 



generation and developing plasma diagnostics techniques. In 
section 2 we describe the instruments and the analysed datasets. 
Section 3 describes the methods used in the analysis. Section 4 
describes the results obtained. Sections 5 and 6 present discus- 
sion and conclusions, respectively. 



2. Instruments and observations 

To study oscillations in the microwave emission from sunspots 
we used the Nobeyama Radioheliograph (NoRH, see Nakajima 
11994b data at 17 GHz. Both correlation curves and interferomet- 
ric images with the spatial pixel size of 4.9 arcsec were analysed. 
The image synthesis and cleaning of compact sources were made 
with the Hanaoka algorithm. The images were synthesised every 
10 s. We use the data obtained in the circular polarisation chan- 
nel, where 3-min oscillations are most evident (e.g. Abramenko 
and Tsvetkov [19851 Shibasaki |200ll Gelfreikh et al. [T999] l. We 
selected 1 1 time intervals during 2005-2010 when only one well- 
developed and isolated sunspot, with a sufficiently strong mag- 
netic field, was present on the solar disk: January 23-25, 2002 
(NOAA 9787), May 02-05, 2005 (NOAA 0756), July 4, 2006 
(NOAA 0898), August 15, 2006 (NOAA 0904) April 29, 2007 
(NOAA 0953) and December 8, 2010 (NOAA 1 1131). Some of 
these events were discussed in our previous studies (Sych and 
Nakariakov|2008l Sych et al. l2009l Sych et al. IMS- 
Figured] shows a typical microwave image, used in the anal- 
ysis. The image of the full solar disk was obtained with NoRH 
at 17 GHz on 3 May 2005 in the circular polarisation chan- 
nel. There is only one active region, NOAA 10756 (marked by 
the rectangle in Fig [T^), with a compact, highly polarised radio 
source of the L-polarity, with the polarisation degree of about 
80% and the brightness temperature of about 127,000 K. Its 
spatial coordinates coincide with the large, symmetric sunspot 
seen in the magnetogram. The strength of the magnetic field 
was estimated with SOHO/MDI. At 01:35 UT on May 3, 2005 
(Fig-IIh) the active region with the centre coordinates (430, -60) 
was mainly unipolar, with prevailing S-polarity (dark areas) and 
a few isolated islands of the opposite polarity (bright areas). The 
value of the magnetic field was about 2,700 G. The polarised 
radio source at 17 GHz (dashed contour lines), corresponding 
to the extraordinary electromagnetic wave, synthesised at 02:45 



UT, basically coincides with the region of the strong magnetic 
field. 

Series of EV and EUV images on May 4, 2005 and 
December 8, 2010 were obtained with the Transition Region and 
Coronal Explorer (TRACE, see Handy and Acton 119991 ) and the 
Atmospheric Imaging Assembly (AIA) of the Solar Dynamic 
Observatory (Title et al. I2006I I. respectively. The TRACE data 
were processed using standard procedures: to remove the back- 
ground and spike-like brightening due to cosmic rays, to cali- 
brate, and where it was necessary to make normalisation of the 
exposure time, and the removal of differential rotation during the 
observation period. The SDO images were obtained from level 1 
data at nine wavelengths: 171, 193,211,304, 131,335,94, 1600 
and 1700A , with the standard processing such as bad-pixel re- 
moval, despiking and flat-fielding and centring the images in a 
temporary data cube. 



3. Methods 

Analysis of transient spatially-varying imaging datasets is a 
nontrivial task that require the use of dedicated tools. In our 
study, we applied the techniques of the wavelet transform with 
the Morlet mother function (Torrence and Compo 119981 ) and 
the Pixelised Wavelet Filtering (PWF, Sych and Nakariakov 
I2OO8I 1. The PWF technique allows us to make narrowband 
maps of power, amplitude and phase distribution of oscillations 
sources at selected harmonics, and obtain their temporal dy- 
namics. The data were prepared with the use the web-based in- 
teractive system for the remote processing of imaging datasets 
(http://pwf.iszf.irk.ru, Sych et al. 120101) . 

Transient increase in the oscillation power, leading to the lo- 
cal increase in the corresponding wavelet coefficients, was anal- 
ysed with the method of wavelet skeletons (Mallat 119991) . The 
method is based on extracting only those coefficients in the 
wavelet decomposition, whose complex amplitudes are high- 
est. Periodicities in a real signal, decomposed to a set of fre- 
quency harmonics of different amplitudes, form local peaks in 
the wavelet spectrum. Mapping their skeletons extracts the most- 
powerful harmonics, and track their time and frequency drifts. 
One can construct two types of skeletons, considering the be- 
haviour of wavelet coefficient extrema either in the frequency 
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Fig. 4. Wavelet analysis of the 3-min oscillation of the microwave signal on May 3, 2005 from 23:55 to 03:15 UT. Panel (a) shows 
the amplitudes obtained by calculating the inverse wavelet transform of the signal in the 120-240 s period range. Panel (b): power 
distribution of the 3-min oscillation trains in the wavelet spectrum. Panel (c): time evolution of the 3-min oscillation power. The 
numbers show the numeration of the trains. Panel (d) shows the wavelet skeleton representing both global (thick lines) and local 
(thin lines) extrema. In the circles we show the signs of the period drifts (positive, negative or without drift). Time is shown in UT, 
power and amplitudes are in arbitrary units. 



or the time domains. In this study, we use the time skeletons, 
applying the following algorithm: 

1 . A power wavelet spectrum with the use of a temporal Morlet 
mother function is constructed. 

2. The frequency and time domains in the wavelet spectrum, 
where oscillations take place, are identified. 



3. All local maxima of wavelet coefficients along the time axis 
in the selected time-frequency region are determined. 

4. The same is repeated for global maxima. 

5. The wavelet skeleton combining the local and global maps is 
obtained. 



4 



Sych et al.: Frequency drifts of 3-minute oscillations in microwave and EUV above sunspots 

Table 1. Parameters of 3-min oscillations: the date, the start and end time of the observation, duration of the observation (in minutes), 
the average period of the spectral power peak, average duration of 3-min oscillation trains, times between the trains, and the fractions 
of the number of trains with positive, negative and variable drifts. 



No 


Date 


Start 


End 


Dtiration 


Main 


Trains 


Time 


Negative 


Positive 


Variable 






UT 


UT 


Observation 


Peak 


Duration 


Separation 


Drift 


Drift 


Drift 






(h,m) 


(h,m) 


(min) 


(min) 


(min) 


(min) 


(%) 


(%) 


(%) 


1 


02.05.2005 


03:05 


06:25 


200 


2.9 


11.7 


13,20,40,61 


53 


34 


13 


2 


03.05.2005 


23:55 


03:17 


200 


2.8 


13.3 


8,12,16,29,52 


60 


32 


8 


3 


04.05.2005 


22:55 


01:45 


170 


2.9 


11.3 


15,30,39,55 


45 


51 


4 


4 


05.05.2005 


00:35 


03:55 


200 


2.6 


10.2 


13,15,21,38,51 


56 


43 


1 


5 


23.01.2002 


23:59 


05:00 


300 


2.8 


10.6 


11,18,30,37,57 


43 


54 


3 


6 


24.01.2002 


01:30 


04:50 


200 


3.0 


11.1 


11,19,27,46,67 


51 


47 


2 


7 


25.01.2002 


23:10 


02:20 


190 


2.8 


11.8 


13,23,32,46 


53 


47 





8 


04.07.2006 


00:45 


05:45 


300 


2.8 


11.6 


16,33,43,58 


43 


51 


6 


9 


15.08.2006 


23:55 


04:55 


300 


2.8 


13.2 


27,36,49,74 


62 


36 


2 


10 


29.04.2007 


23:15 


05:35 


380 


2.8 


11.1 


20,26,42,54 


48 


46 


9 



4. Results 

4.1. Temporal structure of 3-min microwave osciilations 

Figure |2] shows the profile of the NoRH correlation curve ob- 
tained in the circular polarisation channel (R-L) between 22:45, 
3 May 2005 and 06:30 UT, 4 May 2005 with 1-s cadence. The 
relative amplitude of the oscillations is about 3-10%, reach- 
ing in some peaks 17%. The 3-min oscillations, clearly seen in 
the correlation curve, are mainly related to the variation of the 
sunspot radio emission. Also, the correlation curve contains sev- 
eral spikes, corresponding to irregular energy releases. It should 
be noted that the radio flux from the source changed slightly in 
comparison to the full Sun flux, which could be attributed to 
the high level of atmospheric noise. The use of NoRH corre- 
lation curves allows one to magnify the signal associated with 
small-scale sources on the solar disk, and hence to carry out 
preliminary analysis of oscillations without synthesising inter- 
ferometric maps. For the detailed analysis, we chose the time 
interval between 23:55 UT, May 2 and 03:15 UT, May 3, 2005. 
The zoomed part of this interval, shown inside the rectangle in 
Figure |2] clearly demonstrates the presence of well-pronounced 
3-min oscillations. The duration of the interval, 200 min, and the 
cadence time, 1 s, allowed us to study oscillations with periods 
from 3 s to 60 min. 

Spectral wavelet transform of the correlation curve gives us 
its global wavelet spectrum (GWS) shown in Figure[3] The GWS 
spectrum has several well-pronounced peaks, corresponding to 
the periods of about 3 min and also to 8, 12 and 16 min, and 
the long-period components of 30 and 50 min. We would like 
to point out that long-period spectral peaks have already been 
detected in the emission from sunspots, see, e.g. Gelfreikh et al. 
(fT999] l and Chorley et al. (l20T0l l. 

We built a power wavelet spectrum of the correlation curve 
and extracted the periods about 3-min, constructing a narrow- 
band signal and wavelet skeleton shown in Figure |4] It is evi- 
dent that the low-frequency amplitude modulation of the 3-min 
oscillation has a form of repetitive oscillation trains of varying 
amplitude. To obtain parameters of this modulation and visu- 
alise the power distribution of the oscillations on the period-time 
plane, the power wavelet spectrum was made for the narrow- 
band signal filtered in the interval of 120 to 240 s periods. The 
wavelet spectrum is shown in Fig.|4j). Shaded contours show the 
power distribution normalised to the maximum value. The mag- 
nitude contours range from 0.2 to 0.9. Darker areas represent 
more powerful oscillations. Irregular distribution of the power 
of 3-min oscillations in time and periods is clearly seen. The ir- 



regularity in the power is also seen in the main periods of the 
oscillations, asymmetry of the wavelet spectral contours of the 
individual trains, and also in the variations of the train duration 
and time intervals between them. We integrate the wavelet coef- 
ficients over the frequency band, getting the time profile of the 
3-min oscillation power (Fig.|4j;. A significant periodicity in the 
train modulation is clearly seen. In the analysed time interval, 
we find about twelve trains of different maximum power and du- 
ration. The train numbers are indicated on the their tops. 

We would like to point out that the periods evident in GWS 
(Fig-EJ correspond to the characteristic times of the wave trains 
seeing in Fig. |4};. To verify this observation, we found the time 
distance between the beginning and ends of each oscillation 
train. The duration of individual trains range between 7 and 21 
min with the average of about 13 min. These values coincide 
with the evident peaks in the GWS spectrum in the range from 8 
to 16 min. Therefore, it can be asserted that the average duration 
of 3-min trains is about 13 minutes with variation close to 8, 12 
and 16 min. Moreover, the long -period oscillations with the pe- 
riods of 30 and 50 min are possibly related to the characteristic 
time interval of the train occurrence. This possible association 
has also been pointed out in Chorley et al. ( 1201 Oi l. 

The observed amplitude modulation of 3-min oscillations is 
accompanied by the change of the frequency. These changes can 
be visually detected by the behaviour of the power distribution 
in the individual trains in the frequency range (panel b, Fig.|4]i. 
To study this effect, we constructed the wavelet skeleton of the 
NoRH correlation curve (Fig. |4jJ). This allows us to track the 
time and frequency dynamics of both, global and local extrema, 
and to obtain the values of frequency drifts for different phases 
of the sunspot activity. We can see that the period of 3-min com- 
ponent is not stable during the observations. There are repetitive 
frequency drifts with the rate from 3 to 5 mHz/hour, mostly oc- 
curring in the frequency band of 8.3-4.2 mHz, both increasing 
and decreasing the frequency. In the analysed time interval, neg- 
ative drifts prevail, e.g. for trains 2, 3, 7, 8, 10 and 11. There 
are also intervals with positive drifts (trains 5 and 12), and with 
variable drifts (trains 1, 4, 6 and 9). The signs of the drifts are 
indicated in circles in the figure. 

The times of the beginning and ends of the individual trains 
are shown by broken curves in Fig.HJl. Comparing the times of 
the power peaks of the 3-min oscillation trains (Fig. |4]:), and 
their start and end times with the temporal dynamics of the drift 
skeletons, we see that the beginnings of trains and the frequency 
drifts coincide in the majority of cases. The period of 172 s, cor- 
responding to the maximum of the spectral power, determined 
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Time, UT 

Fig. 5. Wavelet skeleton of the NoRH correlation curve in the 90-220 s period band as a function of time, obtained at 17 GHz in the 
polarisation channel on May 04, 2005 03:20 - 05:00 UT. The numbers indicate the times of the narrowband images shown in Fig.|6] 
The mesh strips indicate the prevailing trend of drifts. 



by GWS, usually falls in the middle drift (indicated by the ar- 
rows in Fig. nil). The drift typically begins from shorter periods, 
of about 140 s, and ends at about 220 s. The formation of another 
train, with the frequency significantly different from the previous 
one, is possible before the end of the previous train. Trains 7, 8 
and 9 trains in 02:05-02:20 UT are good examples. In individual 
trains, the linear, systematic frequency drift can be accompanied 
with some deviations, e.g. at 01:47-02:15 UT and 02:27-02:55 
UT 

Similar results were found for other single sunspots anal- 
ysed (02-05.05.2005 NOAA 10756, 23-25.01.2002 NOAA 
9787, 04.07.2006 NOAA 10898, 15.08.2006 NOAA 2006 and 
29.4.2007 NOAA 10953): 3-min oscillations are subject to the 
persistent modulation in all analysed observations. The modula- 
tion depth varies in time by two orders of magnitude. Wavelet 
skeletons show both positive and negative frequency drifts. 
Statistical properties of the parameters of 3-min oscillations de- 
termined by the ten events are shown in Table 1 . 

4.2. Spatial structure of 3-min radio oscillations 

To study the relationship between the frequency drift and spa- 
tial changes in sunspot-associated radio sources, we consider 
sunspot NOAA 10756 during 3:20-05:00 UT on May 4, 2005. 
The GWS of the NoRH correlation curve in the 17 GHz polar- 
isation channel shows well-pronounced 3-min oscillations. The 
wavelet skeleton with the set of global and local extrema of the 
wavelet coefficients is shown in Figure |5] In the analysed time 
interval, there are three well-seen negative frequency drifts, sim- 
ilar to the detected on the previous day (Fig.|4]i. The drift speed 



is about 7 mHz/hour The typical duration of the drift is 40 min. 
The start and end times of consecutive trains are seen to overlap. 
For the times, indicated in Fig. |5] we constructed narrowband 
images of the spatial structure of the 3-min oscillation sources 
with the use of the PWF technique, shown in Figure |6] It is evi- 
dent that there are significant spatial changes in the narrowband 
3-min oscillation sources during the development of the drift: 
new drifts are clearly associated with the appearance of new 
sources of the oscillations. For example the beginning of the drift 
at about 3:40 UT is accompanied by the development of a new 
source labelled by (a) in addition to the already existing source 
(d) that corresponds to the previous act of the frequency drift. 

4.3. Spatial structure of 3-min oscillation trains in EUV 

In order to confirm the results of the frequency-spatial changes in 
3-min oscillations, determined at the transition region in the mi- 
crowave band, we analyse the data obtained at other levels of the 
sunspot atmosphere, in the UV and EUV bands on May 4, 2005 
using TRACE, and on December 8, 2010 using the SDO/AIA. 

4.3.1 . The event on May 04, 2005 observed with TRACE 

EUV images of the active region NOAA 10756 with the cadence 
of about 9 s were obtained with TRACE at 171 A on May 4, 
2005 at 03:30-04:10 UT. After the standard caUbration and pre- 
processing procedure, a data cube of the 50x50 arcsec images of 
the active region was obtained. A movie of the data cube clearly 
shows periodic EUV disturbances propagating along the plasma 
fan from the sunspot outwards. Integrating the data cube over 
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Fig. 6. Time evolution of the spatial sources of 3-min oscillations of the microwave polarisation signal from sunspot NOAA 10756, 
obtained with NoRH at 17 GHz on 2005 May 04, during the 3-min frequency drift 03:20-04:30 UT. The grey colour scale shows 
the power of the oscillation. The dashed curve indicates the umbra-penumbra boundary. The integer numbers label the instants of 
time shown in Fig.|5] The letters label individual sources of the oscillations. 



the spatial coordinates we obtain the time profile of the EUV 
flux from the active region. The wavelet skeleton constructed in 
the 90-220 s period band shows a negative frequency drift of the 
3-min oscillations (Fig.|7] left panel), similar to the drift detected 
in the microwave in the transition region (Fig.|5]l. However, the 
drift speed is about 13 mHz/hour, which is almost two times 
higher than detected at 17 GHz. 

We used the PWF-method investigating the spatial structure 
of sources of 3-min oscillations at diff'erent heights and their re- 
lationship to the frequency drift, at the time about 03:40 UT, 
when significant changes were detected in the microwave band 
(Fig. El right panel). It is evident that the power of 3-min waves 
concentrates in some open plasma structures (which were called 
"wave traces" in Sych et al. 120091 12010b . Although there is al- 
most an order of magnitude difference in the spatial resolution 
of NoRH and TRACE, the 3-min oscillation sources determined 
in EUV have similar spatial location and shape as the sources 
determined in microwaves (see the details labelled by the letters 
a, b, c and d in Fig.|6]and|7] right panel). Like in the microwave 
data, the appearance of new sources of 3-min oscillations coin- 
cides with the beginning of a 30-min train as well as with the 
frequency drift. 

4.3.2. The event on December 08, 2010 observed with 
SDO/AIA 

Active region NOAA 11131 was observed by SDO/AIA on 
December 8, 2010. Sequences of 72x72 arcsec images were 
obtained with the 24-s cadence over the time interval 02:30- 



03:30 UT at 1600, 1700, 171,304, 193,211, 131,335 and 94 A 
wavelengths. This set of observational wavelengths allows us to 
study properties of 3-min oscillations, in particular the frequency 
drifts, from the photospheric up to coronal heights (Figure|8] top 
panel). Applying PWF, we make data cubes of the spatial distri- 
bution of narrowband 3-min oscillations. Snapshots of the spa- 
tial structure of the oscillations are shown in the bottom panels 

of Fig. m 

We see that sources of 3-min oscillations are mainly concen- 
trated in the umbra. For lower layers of solar atmosphere, i.e. at 
1700 A, the sources have the shape of symmetric patches that 
fill almost uniformly the central part of the umbra. In the chro- 
mospheric level, i.e. at 304 A, we see the fragmentation of the 
oscillation source that gains a spiral shape, filling up the whole 
umbra. In the coronal level (171 A), further deformation of the 
oscillation source takes place: a part of the source has a radial 
sector structure, while the other part is stretched along coronal 
plasma loops that form the active region fan. A similar shape 
of the oscillation source is seen in higher-temperature emission 
lines, but the power of 3-min oscillations in this emission is sig- 
nificantly lower. 

Information about frequency drifts and wave propagation 
can be obtained from time-distance plots, constructed for the 
time variation of the emission intensity along a certain one- 
dimensional spatial slit. For example, such a plot for a horizontal 
slit passing through the centre of the sunspot in the 304A data 
cube is shown in Figure [TO] The wave structure of 3-min oscil- 
lations, propagating symmetrically from the sunspot centre, the 
consequent decrease in the propagation speed and damping at 
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Fig. 7. Left panel: Wavelet skeleton of 3-min oscillations observed in the EUV emission over sunspot NOAA 10756 on May 04, 
2005 at 03:30-04:08 UT. The negative trend of the frequency drift is shown by the mesh strip. Thick lines correspond to the global 
maxima of wavelet coefficients, and the thin lines to the local ones. Right panel: Overlap of the narrowband image of 3-min EUV 
oscillation sources (shown by white contours) on the integral EUV image at 03:40 UT. The letters a, b, c and d correspond to the 
fine 3-min oscillation radio sources in Fig.|6]at 03:45 UT. The arrows indicate the apparent direction of the wave propagation. 



the umbra-penumbra boundary, are clearly seen. Analysis of the 
time characteristics of the signals was performed by the same 
fashion as of the TRACE data (Sec. 14.3. Il l separately for each 
observational wavelength. 

Figure |9] shows the amplitude and phase modulation of 3- 
min oscillations for three wavelengths, 1700, 304 and 171 A , 
which correspond to three different heights, and their wavelet 
skeletons. As in our previous findings, the oscillation trains last 
typically from 8 to 20 minutes (with the average of about 12 
min), showing partial overlap. Times between individual trains 
are about 30 minutes. Wavelet skeletons show that the rising 
phase of the trains is accompanied with negative frequency 
trains. The average period of 3-min oscillations is 155 s at 1700 
A, 153 s at 304 A and 162 sec at 171 A . Negative frequency 
drifts dominate in all observational channels. The drift speed at 
the photospheric level (1600 A, 1700 A) is about 4-5 mHz/hour, 
at the chromospheric level (304 A) about 5-7 mHz/hour, and at 
microwaves (NoRH, 17GHz) 5.6-8 mHz/hour and in the corona 
(171, 193, 211 and 131 A) about 11-13 mHz/hour. Thus, the 
absolute value of the drift speed increases by almost two times 
during the propagating from the sub-photospheric levels to the 
corona. In higher temperature emission (335 and 94 A) the 3- 
min signal is weak, and hence it is impossible to get reliable 
information about the drifts. 



5. Discussion 

Using high-quality EUV and microwave data obtained with 
SDO/AIA, TRACE and NoRH and modern analytical tech- 
niques (PWF-method and wavelet skeletons), we studied com- 
plex temporal dynamics of 3-min oscillations above sunspots. 
The dynamics includes variations in height, spatial, temporal 
and frequency structure of the oscillations. Spatially-integrated 
signals obtained at different heights and with different instru- 
ments show common behaviour: the oscillations experience am- 



plitude and frequency modulation, which is accompanied with 
the change of the spatial sources. The amplitude modulation has 
the form of wave trains, of the duration from 7 and 21 min and 
the repetition period of about 30-60 min. The frequency modu- 
lation has the form of frequency (or period) drifts in the 2-4 min 
period range. Typically, the beginning and end times of the drifts 
coincide with the trains of the oscillations, which indicates their 
phenomenological relationship. Mainly negative drifts (the grad- 
ual increase in the period of 3-min oscillations) are observed. 
Sometimes two drifts co-exist, i.e. during the end of the pre- 
vious drift, when the frequency is in the low-frequency part of 
the 3-min spectrum, near 200 s, a new drift appears near 160 s. 
Development of new 3-min oscillation trains is typically accom- 
panied with the appearance of new spatial sources (the regions of 
the enhanced power of 3-min oscillations in the horizontal plane) 
in the sunspot atmosphere. In the corona, spatial sources of 3- 
min oscillations develop along certain magnetic flux tubes of 
the magnetic fans above the sunspot. Narrowband time-distance 
plots show that in those structures 3-min oscillations have a form 
of outwardly propagating waves. This behaviour is well consis- 
tent with recent numerical modelling (Botha et al. 1201 11 1 and 
the analysis of imaging and spectral EUV data (e.g. De Moortel 
120091 Verwichte et al. 1201 Oi l. Also, there are coronal sources of 
another type: compact, of a symmetric, patch-like shape, which 
can perhaps be associated with the magnetic flux tubes with a 
narrow angle to the line-of-sight. 



5.1. A possible mechanism for frequency drifts 

One of the possible mechanisms explaining the observed prop- 
erties of 3-min oscillations oscillation is the response of the 
sharply stratified atmosphere to pulses coming from lower re- 
gions. The pulses can be excited by magnetic reconnection, 
granulation or other explosive or impulsive events taking place in 
sub-photospheric, photospheric or lower chromospheric regions. 
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Fig. 8. Top panels: Images of NOAA 11131 at 1700 A (photosphere), 304 A (chromosphere) and 171 A (corona) on December 8, 
2010 at 02:30 UT. Bottom panels: spatial distribution of 3-min oscillations, obtained by filtering in the 5.2-9.5 mHz band. Thin 
dotted lines indicate the umbra-penumbra boundary at the photospheric level. The thick dashed line in the narrowband image at 304 
A shows the slit of the time-distance plot. 

Solar Dynamic Observatory (SDO), December 08, 2010 
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Fig. 9. Top panels: Time profiles of 3-min oscillation trains in the emission intensity of the 1700, 304 and 171 A signals obtained 
with SDO/AIA on December 8, 2010 at 2:30-03:30 UT. Envelopes of the trains are shown by the thick curves. Bottom panels: 
wavelet skeletons of the signals. The continuous line shows the period corresponding to the maximum spectral power. 
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Fig. 10. The time-distance plot of the active region NOAA 11131 
on December 8, 2010 obtained by SDO/AIA at 304 A . The in- 
tensity oscillations are presented in the logarithmic scale. The 
white contours show the narrowband wave fronts in the 5.2- 
9.3 mHz band. The spatial slit is positioned across the sunspot, 
through its centre, and is shown in Fig. [8] 



Consider a gravitationally stratified plane atmosphere with a 
vertical magnetic field. We assume that po{z) and po(z) are the 
plasma pressure and density, respectively; and the temperature. 
To, is constant. Then sound speed, Cj - ^JypQ/pQ, is also con- 
stant. The pressure scale height is 



H = 



fhg ' 



(1) 



where kg is the Boltzmann constant, g is the gravitational accel- 
eration and m is the mean molecular weight. This model of the 
solar lower atmosphere is obviously over-simphfied, but we use 
it for illustration purposes only. 

Propagation of linear sound waves in the stratified medium 
is governed by the Klein-Gordon equation, which implies the 
acoustic cut-off frequency 



Uc - 



2H 



(2) 



Waves with frequencies higher than cjc can propagate upwards, 
while the lower-frequency waves are evanescent. However, evo- 
lution of an initially broadband pulse leads to a formation of an 
oscillatory wake (Lamb |T908l Rae & Roberts [19821 ): 



Ml — UqJo 



2^ 

Z 

4^2 



where ui is the vertical component of the velocity, uq is the initial 
amplitude, Jq is the Bessel function of the zeroth order, and 'H 
is the Heaviside function. The time is measured from the time 
when the pulse was launched. 

As the chromospheric cut-off period is about 3-min then 
the observed oscillations can be explained as wakes behind the 
pulses propagating along the vertical magnetic field (Fleck and 
SchmitzJ991 , Kalkofen et al. 1 19941 Sutmann and Ulmschneider 
119951 Kuridze et al. I2009I I. However, each active region may 
have a fine structure in the horizontal direction, in the form of 
flux tubes with different plasma temperatures, as the strong mag- 
netic field prohibits thermal conduction across the field. Then, 
the pulses propagating in different parts of the active region may 
set up the wakes with different oscillating periods. The radio ob- 
servations may not resolve the fine horizontal structure of the 



oscillations due to insufficient spatial resolution. Then the oscil- 
lations with different frequencies will be superimposed, and the 
resulted observed oscillation may have a complex structure. 

In order to illustrate this scenario, in Figure [TT] we plot the 
linear superposition of two wakes, which could be formed by 
two acoustic pulses propagating along two different flux tubes 
in the umbra. The ratio of the cut-off frequencies in these two 
regions is taken to be 0.8. The second pulse was launched 810 
s later than the first pulse. We also show the wavelet spectrum 
of the signal, time variation of the oscillation power, and the 
wavelet skeleton with the global maximum. The time interval 
between the wave trains is about 14 min, which actually is the 
interval between the consecutive pulses. The oscillation period 
slightly increases (decreases) before (after) the maximum of the 
oscillation amplitude. Thus, the behaviour of the modelled oscil- 
lation is very similar to the observed behaviour of 3-min oscilla- 
tion trains. 



6. Conclusions 

We summarise results of this study as follows: 



1. 



We analysed 3-min oscillations in the NoRH polarisation 
correlation curves at 17 GHz, obtained for eleven time in- 
tervals when there was only one powerful sunspot on the 
solar disk. The relative amplitude of 3-min oscillations was 
about 3-8% reaching 17%. Spectral analysis demonstrated 
the presence of significant peaks in the range of 2-4 min with 
a power concentration near 3 min period. The 3-min oscilla- 
tions experience amplitude modulation in the form of peri- 
odic wave trains, lasting for 8-20 min with 13 min in average. 
The repetition rate of the trains is 20-40 min with 30 min in 
average. There were also low-frequency peaks in the spec- 
trum, corresponding to the periods of 10-20 min, and 30-60 
min. Apparently these long-period oscillations coincide with 
the characteristic times of the 3-min oscillation trains. 
In addition, 3-min oscillations are subject to frequency mod- 
ulation: their period varies in the range 90-240 s. During the 
oscillation train, the frequency modulation has often a form 
of a steady frequency drift. The beginning and end of the 
drift coincide with the beginning and end of the amplitude 
modulation. In the majority of cases, the negative drift, when 
the 3-min oscillation period gradually increases, is seen. In 
different trains of 3-min oscillations observed in the same 
sunspot, there can be both positive and negative drifts. The 
start and end times of consecutive drifts can overlap. This be- 
haviour is seen at all levels of the sunspot atmosphere. The 
speed of the drift is 4-5 mHz/hour in the photosphere, 5-8 
mHz/hour in the chromosphere, and 11-13 mHz/hour in the 
corona. Thus, apparently, the drift speed grows with height. 
In addition to the gradual drift, the instant frequency of 3-min 
oscillations experience positive and negative fluctuations. 
Analysis of the spatial structure of 3-min sunspot oscilla- 
tions observed in microwaves and EUV shows that the start 
of drifts coincides with the appearance of sources of the os- 
cillations (the regions with the enhanced power of 3-min os- 
cillations). 

Observed spatial-frequency-temporal properties of 3-min os- 
cillations can be interpreted in terms of two simultaneously- 
operating mechanisms. The first effect is connected to the 
negative frequency drift due to the dispersive evolution of 
upward propagating sub-photospheric pulses, guided by the 
sunspot magnetic field. The repetition rate of the pulses (30- 
60 min) is similar to the long-period oscillations of sunspots. 
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Fig. 11. A model of the observed behaviour of 3-min oscillations: superposition of two oscillatory wakes formed in two different 
magnetic flux tubes in a sunspot's umbra. Panel (a): the signal; panel (b); its wavelet spectrum; panel (c): time dependence of the 
oscillation power; panel (d): wavelet skeleton of the signal. 



which suggest their possible association. The second effect 
is the spatial separation of the propagating pulse by several 
individual magnetic flux tubes, of different curvature, tem- 
perature and line-of-sight angles. In particular, the frequency 
drifts may take place due to the changes in the tube curvature 
with height. 

Acknowledgements. The work was supported by the grants RFBR 10-02-00153 
and 1 1-02-10000-k. RS is supported by the Chinese Academy of Sciences 
Visiting Professorship for Senior Intemational Scientists No. 2010T2J24. TZ 
is supported by the Austrian Fonds zur Forderung der wissenschaftlichen 
Forschung (project P21 197-N16) and the Georgian National Science Foundation 
(grant GNSF/ST09/4-310). The TRACE, SDO and NoRH data were obtained 
from the TRACE, SDO and NoRH databases. The authors acknowledge the 
TRACE, SDO and NoRH consortia for operating the instruinents and perforin- 
ing the basic data reduction, and especially for the open data policy. Wavelet 
software was provided by C. Torrence and G. Compo, and is available at URL: 
[http://paos. Colorado. edu/research/wavelets/| 



References 

Abramenko, V. I., Tsvetkov, L. I., 1985, Bull. Crim. Astrophys. Obs. 73, 53 
Balthasar, H., Martinez Fillet, V., Schleicher, H., 1998, Solar Phys. 182, 65 
Bard, S., Carlsson. M., 1997, SOHO Workshop: The Corona and Solar Wind 

Near Minimum Activity ESA SP-404, 189 
Bellot Rubio, L. R., Collados, M., Ruiz Cobo, B., Rodriguez Hidalgo, I., 2000, 

ApJ, 534, 989 
Bogdan, T. J., 2000, Sol. Phys., 192, 373 

Botha, G. J. J., Arber, T. D., Nakariakov, V. M., & Zhugzhda, Y. D. 201 1, ApJ, 
728. 84 

Brkovich, A., Ruedi, I., Solanki, S. K. et al., 2000, A&A, 353, 1083 

Chorley N., Hnat, B., Nakariakov, V. M., Inglis, A. R., & Bakunina, I. A. 2010, 

A&A, 513, A27 
De Moortel, 1. 2009, Space Sci. Rev., 149, 65 
Elliott, 1., 1969, Solar Phys., 6, 28 



Fleck, B., Schmitz, F., 1991, A&A, 250, 235. 
Fludra. A.. A&A, 2001, 368, 639 

Gelfreikh G. B., Grechnev, V., Kosugi, T. et al., 1999, Sol. Phys. 185, 177 
Gelfreikh, G. B.; Shibasaki, K., Nagovitsyna, E. Yu. et al., 2004, lAU 

Symposium, 223, Cambridge University Press, 245 
Gurman, J. B., Leibacher, J. W., Shine, R. A. et al., W., 1982, ApJ., 253, 939 
Handy B. N., Acton, L. W., 1999, Solar Physics, 187, 229 
Harrison, R. A.. Lang, J., Brooks, D. H. et al., 1999, A&A, 351, 1115 
Kalkofen, W., Rossi, P, Bodo, G.et al., 1994, A&A, 284, 976 
Kobanov, N. I., Kolobov, D. Y, Chupin, S. A., & Nakariakov, V. M. 201 1, A&A, 

525. A41 

Kobanov, N. I.. Kolobov. D. Y, & Chupin, S. A. 2008, Astronomy Letters, 34, 
133 

Kuridze, D., Zaqarashvili, T. V., Shergelashvili, B. M. et al., 2009, A&A, 505, 
763 

Lamb, H., 1908, Proc. Lond. Math. Soc, Ser, 2, 7, 122 
Lites, B. W., 1986, ApJ, 301, 1005 

Lites. B. W., 1992, in Sunspots: Theory and Observations, ed. J. H. Thomas N. 

O.Weiss (NATO ASJ Ser C, 375; Dordrecht: Kluwer), 261 
Leibacher J. W., Stein, R. F., 1981, In The Sun as a Star, Jordan, S. (ed.), NASA 

SP-450. 263 

Leibacher J., Gouttebroze, P, Stein, R. F., 1982, ApJ, 258, 393 
Liu, S.Y, 1974. ApJ, 189, 359 

Mallat, S., 1999, A Wavelet Tour of Signal Processing, 2nd edn., San Diego 

Acadeinic Press 
Murawski, K., Zaqarashvili, T. V., 2010, A&A, 519, A8 
Muzy J. F. Bacry, E. Arneodo, A. 1993, Physical Review E, 47. 875 
Nakajima, H., Nishio, M., Enome, S., et al. 1994, IEEE Proceedings, 82, 705 
Nakariakov, V. M., Inglis, A. R., Zimovets, 1. V. et al., 2010, Plasma Physics and 

Controlled Fusion, 52, 12 
Nindos, A., Alissandrakis, C. E., Gelfreikh, G. B. et al., 2002, A&A, 386, 658 
Norton, A. A., Ulrich, R. K., Bush, R. I. et al., 1999, ApJ, 518, L123 
Orrall, RQ., 1966, ApJ, 143, 917 

Percival, D.B. Walden A.T., 2000, "Wavelet methods for time series analysis", 

Cambridge University Press 
Rae I. C, Roberts, B., 1982, ApJ, 256, 761 



11 



Sych et al.: Frequency drifts of 3-minute oscillations in microwave and EUV above sunspots 

Rendtel, J., Staude, J., Curdt, W., 2003, A&A, 410, 315-321 

Rioul, O, Vetterli, M., 1991, "Wavelet and Signal Processing", IEEE Signal 

Processing Magazine 
Shibasaki, K. 2001, ApJ, 550, 1113 
Sutmann, G., Ulmschneider, P, 1995, A&A, 294, 232 
Sych, R. A., Nakariakov, V. M., 2008, Solar Physics, 248, 395 
Sych, R., Nakariakov, V. M., Karlicky, M.et al., 2009, A&A, 505,791 
Sych, R. A., Nakariakov, V. M., Anfinogentov et al., 2010, Solar Physics, 266, 

349 

Title, A. M., Hoeksema, 1. T. Schrijver, C. J., 2006, 36th COSPAR Scientific 

Assembly, CD ROM 2600 
ToiTence, C, Compo, G.P, 1998, Bull. Am. Meteorol. Soc. 79, 61 
Tziotziou, K., Tsiropoula. G.. Mein. P, 2002, A&A, 381, 279 
Verwichte, E., Marsh, M., Foullon, C, et al. 2010, ApJ, 724, L194 
Zaqarashvih, T. V., Murawski, K., Khodachenko, M. L., Lee, D., 2011, A&A, 

529, A85 

Zhugzhda, Y. D., Balthasar, H., Staude, J., 2000, A&A, 355, 347 



12 



